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Abstract 

The dynamics of a stabihzed radion in the Randah-Sundrum model (RS) with two 
branes is investigated, and the effects of the radion on electroweak precision observ- 
ables are evaluated. The radius is assumed to be stabilized using a bulk scalar field as 
suggested by Goldberger and Wise. First the mass and the wavefunction of the radion 
is determined including the backreaction of the bulk stabilization field on the metric, 
giving a typical radion mass of order the weak scale. This is demonstrated by a per- 
turbative computation of the radion wavefunction. A consequence of the background 
configuration for the scalar field is that after including the backreaction the Kaluza- 
Klein (KK) states of the bulk scalars couple directly to the Standard Model fields on the 
TeV brane. Some cosmological implications are discussed, and in particular it is found 
that the shift in the radion at late times is in agreement with the four-dimensional 
effective theory result. The effect of the radion on the oblique parameters is evaluated 
using an effective theory approach. In the absence of a curvature-scalar Higgs mixing 
operator, these corrections are small and give a negative contribution to S. In the 
presence of such a mixing operator, however, the corrections can be sizable due to the 
modified Higgs and radion couplings. 
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1 Introduction 



Extra dimensional theories where standard model fields are localized on a brane ^ ^ 
S B 0' 01 have recently attracted a lot of attention, since such models have several distinct 
features from ordinary Kaluza-Klein (KK) theories. In particular, Randall and Sundrum 
@] presented a simple model based on two branes and a single extra dimension, where the 
hierarchy problem could be solved due to the exponentially changing metric along the extra 
dimension. In order to obtain a phenomenologically acceptable model, the radion field (which 
corresponds to fiuctuations in the distance of the two branes) has to get a mass, otherwise 
it would violate the equivalence principle 0, and also result in unconventional cosmological 
expansion equations ||T^, |1T|. The simplest mechanism for radius stabilization has been 
suggested by Goldberger and Wise |jl2[, who employed an additional bulk scalar which has 
a bulk mass term and also couples to both branes (for another issues related to the radion 



potential see |]T3|)- Both the cosmology and collider phenomenology crucially depend on the 
mass and couplings of the radion. In particular, there is no radion moduli problem if the 
radion mass is O(TeV) and its couplings to Standard Model (SM) fields is 0(TeV-^). This 
is also the most favorable scenario for discovering the radion at a future collider. 

In fact it was shown in p^ , |T^, that the radion will have the above properties for the 



Goldberger- Wise scenario. However, the calculation of [|1J, [T^ were using a naive ansatz 
for the radion field which ignores both the radion wavefunction and the backreaction of the 
stabilizing scalar field on the metric. The validity of this approximation has recently been 
questioned 



Therefore in this paper we analyze the coupled radion-scalar system in detail from the 
5D point of view. We derive the coupled differential equations governing the dynamics of the 
system, and find the mass eigenvalues for some limiting cases. Due to the coupling between 
the radion and the bulk scalar, we find that there will be a single KK tower describing the 
system, with the metric perturbations non-vanishing for every KK mode. This implies that 
the Standard Model fields localized on the TeV brane will couple to every KK mode from 
the bulk scalar, and this could provide a means to directly probe the stabilizing physics. 

Using the coupled equations for the radion-scalar system, we analyze the late-time be- 
havior of the radion in an expanding universe, and find that the troubling 55-component of 
Einstein's equation just determines the shift of the radion. This shift completely agrees with 
the shift obtained in using the 4D effective theory. 



Given that we have established that the radion mass is O(TeV) and that its couplings 
to SM particles is 0(TeV~^), it is reasonable to consider its effects on SM phenomenology. 
Some direct collider signatures for the radion and loop corrections have been discussed in 



l|, ^ In the second half of this paper the effects of the radion on the 



oblique parameter^ are calculated using an effective theory approach similar to . Since 
in the RS model the radion is the only new state well below the TeV scale, a low-energy 
effective theory including only the radion and SM fields is used. The effects of other heavy 
modes are accounted for by including non-renormalizable operators at the cutoff scale. In 



*Loop effects for tlieories witli large extra dimensions have been analyzed in 
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the absence of a curvature-scalar Higgs mixing operator, the corrections from the radion are 
small, but give a negative contribution to S. In the presence of such a mixing operator the 
corrections could be sizable due to the modified radion and Higgs couplings. 

This paper is organized as follows: in Section ^ we review the Randall- Sundrum model 
and radius stabilization by bulk scalar fields. We also summarize the explicit example of 
p6| which we will be using for our explicit computation of the radion mass and couplings to 
SM fields. In Section ^ we present our ansatz for the coupled metric and scalar fluctuations 
based on the analysis of |27, ^ of the radion without a stabilizing potential. We will 
derive a single ordinary differential equation, whose eigenmodes will yield the KK modes for 
the radion-scalar system. In Section |^ we analyze the generic properties of this equation. 
In the general case we find that the system is not described by a hermitian Schrodinger 
operator. However, we identify a convenient limit, in which the differential operator is in 
fact hermitian, and the eigenf unctions are manifestly orthogonal. In Section ^ we analyze 
the eigenf unctions in this limit, and find the approximate masses for the KK tower. In this 
analysis, the back reaction of the metric is neglected, which results in the lightest mode still 
being massless. The effect of the back reaction on the lightest mode is taken into account in 
Section ^ where we find the mass of the radion to be of the order (but slightly lighter) than 
the weak scale. In Section |^ we discuss the couplings of the radion and the KK tower to 
SM fields on the brane. We find that the radion coupling is exactly agrees with the results 
in [ll4| , [T5| , while the couplings of the other KK modes of the scalar field are suppressed by 
the mass of the given mode, and is proportional to the backreaction of the metric due to the 
scalar background. In Section ^ we demonstrate that in an expanding universe the shift in 



the radion at late times agrees with the 4D effective theory result obtained in WM. Having 



established the mass and coupling of the radion, we write an effective Lagrangian in Section 
^ without any specific mechanism of radius stabilization and neglecting the contributions 
of the KK modes. In Section |1^ we add a curvature-scalar Higgs coupling to the effective 
Lagrangian, and discuss how the couplings are modified. Then, in Section |ll| we calculate 
the Feynman rules in a general gauge. These allow us to compute the oblique parameters via 
one-loop vacuum polarization diagrams with radions in the loop in Section |12|. The radion 
correction is log divergent (unlike the Higgs), and so we also write the nonrenormalizable 
operators at the cutoff scale that provide the necessary counterterms. The size of the new 
contributions are shown for various cases in several figures in Section |12.1| . We also estimate 
limits on the radion mass as a function of the cutoff scale in Section Finally, we conclude 
in Section |14. 



2 Review of the Randall- Sundrum Model and the 
Goldberger-Wise Mechanism 

Randall and Sundrum presented a very interesting proposal for solving the hierarchy problem 
0]. By introducing a fifth dimension where the bulk geometry is anti-de Sitter, a large 
hierarchy between the Planck scale and the TeV scale is obtained with only a mild finetuning. 
Two branes are introduced, located at the boundaries of the anti-de Sitter space. By tuning 



2 



the bulk cosmological constant A = — 6A;^//t^, the tensions Vp and Vr on the Planck and 
TeV branes, respectively, such that Vp = —Vp = 6k /k^ (where is the 5D Newton constant 
related to the 5D Planck mass by = 1/2M'^) one obtains a 4-D Poincare invariant solution. 
The metric is then 

^^2 _ g-2fcy^^^^^M^^i. _ ^^2^ (2.1) 

where the Planck brane and TeV branes are located at y = and y = tq. For a moderate 
choice of kr^ ~ 0(50), a large hierarchy between the Planck scale and the weak scale is 
generated. 

Since this solution is obtained for any value of tq, some mechanism is required to fix 
ro ~ 50/A; as opposed to some other value of tq. This must also be done without introducing 
any large finetuning. Further, small shifts in the separation between the two branes do 
not change the energy, and so are described in an effective theory by the fluctuations of a 
massless particle, the "radion" . This particle couples like a Brans-Dicke scalar and must be 
massive to recover ordinary 4-D Einstein gravity |14| 

One way to achieve these requirements is to introduce a bulk scalar field that has a 
bulk potential V{(j)) [1^. To stabilize the brane distance, potentials Xp,t{4>) on the Planck 



and TeV branes respectively are also included. The competition between the brane and 
bulk Lagrangians generates a vacuum expectation value (vev) for 0, which results in a 4-D 
vacuum energy that depends on tq. For a simple choice of polynomial potentials a large 
hierarchy is then easily obtained with a mild finetuning [1^, and the resulting mass for the 
radion is O(TeV) IH, IH. 

The phenomenology of the radion depends on the strength of its coupling to the brane 
fields. Using the following naive ansatz to describe the radion h{x) 

ds^ = e-2^l^l''(")ds2 - h{xfdy^, (2.2) 

1^ and |15[ computed the normalization of the radion kinetic term to be 



3 _9z,«„ kv. 
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-e --^-^{dhf. (2.3) 
4 

Fields living on the TeV brane couple to the radion through the induced metric, with an 
interaction 

^6(^)TrT,. = ^TrT,. (2.4) 

where r is the canonically normalized radion, Avk = Mpie~^'^° ~ O(TeV), Mpi = (1 — 
g-2fcro^^^^^2'j ^ 1/(A;k^), and T^,y is the physical energy-momentum tensor of the TeV brane 
fields. It is then clear that the radion couples as ~ 1/TeV to the Standard Model fields. 
Obtaining an acceptable phenomenology then requires that the radion mass is O(TeV), which 
is easily satisfied by the Goldberger-Wise mechanism. 

The phenomenology of the radion is then crucially dependent on the normalization of 
the kinetic term. In fact, in the computation leading to the O(TeV)^ prefactor in (|2.3| ) 
there is a cancellation between two terms of 0{Mpi). The origin of this cancellation remains 
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somewhat mysterious, and the absence of this cancellation would clearly lead to different 



predictions. In [0 it was pointed out that there are additional contributions to the radion 



kinetic term not included in 111, ffHJ. In particular, the profile of the stabilizing field depends 



on ro, and so a small change in tq — > b{x) distorts the background field. It was found that 
this results in an 0{Mpi) correction to the radion kinetic term, thereby drastically changing 
the phenomenology of the radion. 

We review the resolution of this issue in the first part of this paper. Some of the results 
presented in Sections ^|-^ are already contained in the work by Tanaka and Montes [p9 | 



even though the results of this paper were obtained independently of Ref . p9| . We explicitly 



determine the wavefunction of the radion when there is a stabilizing mechanism. We find 
that the radion mass is typically O(TeV). In the limit that the backreaction of the stabilizing 
fields on the metric is small, we find that the correction of the stabilizing field to the radion 
kinetic term is subdominant to the gravitational contribution. We also find that once the 
stabilizing field has a non-zero vev, the Kaluza-Klein (KK) tower couples directly to the 
brane world fields, with 1/TeV normalization, and amplitude depending on the size of the 
backreaction. 

The action we consider is |^ 

-M'l d'>x^n + 1 d''x^{^V<pV<P-V{<P)^- j rf^Xv^Ap(0)- j d'x^.Xri^), (2.5) 

where is the induced metric on the branes. The background vev for (p and background 
metric that preserve 4-D Lorentz invariance is 

0(x,y) = 0o(z/), (2.6) 
ds^ = e-^'^r]^^dx^'dx'' - dy^. (2.7) 
The Einstein equations are then 

Rab = n^fab = (Tab " ^^7a6^/'''T,rf) , (2.8) 

with = 1/ (2M^). For this background the scalar and metric field equations are 

4A''-A" = -^K(0o)-yE^^(</'o)%-l/^), (2.9) 

A" = ^-^nM (2.10) 
4 = 4^Vi + ^ + i:^%-!/.). (2.11) 



Here primes denote d/dy, and we reserve to denote derivative with respect to the comov- 
ing 4-D spacetime coordinates x'^. The boundary equations for A and (po are obtained by 
matching the singular terms in the above equations. This gives 

[A% = ^\.{<Po), (2.12) 



*The action is integrated over the circle rather than the hne segment. 
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(2.13) 



For analytical solutions we use an approach presented in j^G], |30|. A particular class of 
potentials V is considered which can be written in the form 



(2.14) 



Then a solution to the following first order equations 



IdW 



A' 



K 

"6 



W^(0o) 



(2.15) 



automatically solves both the Einstein and scalar field equations, once the appropriate 
boundary conditions are solved. The virtue of this method is that for simple choices of 
W it is possible to also solve for the backreaction of (po on the metric. This will be important 
for us, since we find that only after including the backreaction of the stabilizing field does 
one find that the radion acquires a mass. 

In particular, to obtain some analytic results the following superpotential will be 
used: 



W(0) 



6k 



— uc 



with brane potentials 

A(0)± = ±W{<P±) ± W\<P±){<P - 0±) + 7i( 



(2.16) 



(2.17) 



Here +/— refer to Planck/TeV brane. The solution is ||26 



-uy 



My) 



^ 12 



The separation distance tq is then fixed by matching 
gives utq = ln0p/0T. So the quantity 



(2.18) 
(2.19) 

at and tq to (pp and 0t which 

(2.20) 



is not a (hierarchically) small number, since both (pp and 0t are 0{Mp{'^). This combination 
will appear later in the expression for the radion mass. Also for future reference, since 
the backreaction corresponds to the second term in A, the limit of a small backreaction is 



<^ 1, and M > 0, but with (pp/cpT =constant, so that u is kept constant. 
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3 The Coupled Field Equations 



When 00 = there is always a static solution independent of the value of tq.* The small 
fluctuations in the relative position between the two branes then describe a massless particle 
("the radion"), and its wavefunction is G{x,y) = 2F{x,y) = 2ke'^^yR{x) and where 
nR = 0. Since the coupling of the radion to the Standard Model fields is ~ 1/TeV |T5| , 
obtaining an acceptable phenomenology requires that this radion acquires a mass. 

We therefore consider the spectrum of perturbations about the above background which 
stabilizes the inter-brane separation. A general ansatz to describe the spin-0 fluctuations is 

(j){x,y) = (f)o{y) + ^{x^y) (3-1) 
ds^ = e-^^-^^^''^y\^,dx^'dx^ - (1 + y)fdy\ (3.2) 

In order to describe all gravitational excitations of the model, one would need to add also 
the degrees of freedom in the graviton, by replacing rj^y rj^y + /i^J. One can show that the 
Einstein equations with this replacement will have the radion and the graviton decoupled. 
This metric ansatz (|3.2|) (together with the two equations ( p. 12 ) and G = 2F which we 



will shortly derive) fixes our gauge choice. One can show, that the effect of the remaining 
gauge transformations that preserve the form of (|3.2| ) and ( p.l2| ) just amount to a 4D gauge 



transformation on the graviton field h^,^ and can be used to impose a convenient 4D gauge 
for the graviton. In the following we will only concentrate on the radion field. Using this 
ansatz the Einstein and scalar field equations are linearized to obtain some coupled equations 
for F, G and if. The linearized Einstein equations are 

SRab = K'^Sfab- (3.3) 

Inspecting the 6Rfj,iy equation one immediately concludes that G = 2F. For 

6R^, = ■■■ + 2d^d,F - d^d^G + • • • (3.4) 

where the ellipses all contain terms ~ rj^^. Since to linear order in the perturbations the 
sources 5T^y are also all ~ r^^j,, the 9^(9y term in SR^^ term must vanish. This gives G = 
2F + c. However in the limit F — 0, or G — we should recover the background solution, 
so c = 0. In what follows we set G = 2F. Then the coupled field equations are 

SR^^u = Vt^unF + e~^%,[-F" + 10A'F' + 6A"F -24:A'^f) , (3.5) 
5/?m5 = Sd^F' -QA'd^F, (3.6) 
6R55 = 2e'^^nF + 4F" -IGA'F'. (3.7) 

The source terms are 

Sf,, = ~e-^%AV'{My^-2V{(Po)F) 



* After two finctunings which arc independent of tq. But only one finetune remains after radius stabiliza- 

tion p,ra. 
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-^e"^ V E (^^^ - 4A.(0o)f) 6{y - y,), (3.8) 
5f^5 = (l^od.V, (3.9) 




if + 2K{<Po)F]5iy-y.:). (3.10) 



The linearized scalar field equation is 

-6Kf"-4|^F. (3.11) 

Notice that the i?^5 may be integrated immediately to obtain 

(l)',ip=^{F'-2A'F). (3.12) 

An integration constant k{y) has been sent to zero since we require that the fluctuations F 
and If are also localized in x. This equation (p.l2|) together with the metric ansatz ( |3.2| ) 



fixes our gauge choice. One can show, that the effect of the remaining gauge transformations 
that preserve the form of ( p.2| ) and ( |3.12| ) just amount to a 4D gauge transformation on the 
graviton field /i^j/ and can be used to impose a convenient 4D gauge for the graviton. 

These equations must be supplemented by the boundary conditions for F and cp on the 
two branes. These are obtained by identifying the singular terms in above equations. A 
priori the Einstein equations give two boundary conditions for each wall. It is however 
straightforward to show that one of them is trivially satisfied once A satisfies the jump 
equation (|2.12|). The two remaining boundary equations are 



[F'] = ^H<Po)F + ^^{Mv. (3.13) 

[^% = ^^{<Po)^ + 2^F (3.14) 

Upon using the jump equations for the background the first equation is seen to be equivalent 
to ( |3.12| ) and so provides no new constraints. Then only the second boundary condition 
must be implemented. A convenient limit will at times be considered in this paper. The 
second boundary condition simplifies in the limit of a stiff boundary potential. Namely, if 
d'^Xi/dcf)'^ ^ 1 then the second boundary condition is just if\i = 0. Then in this limit the 
first boundary condition is just 

{F' -2A'F)\i = 0. (3.15) 
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A single equation for F is obtained as follows. One considers the combination e^^SR^u + 
5i?55 in the bulk. The point of this combination is to eliminate terms of the form V{(j)Q)!f. 
This leaves a bulk equation for F and if' only: 



e^^aF + F" - 2A'F' = 

3 



(3.16) 



One then eliminates if' in favor of F using ( p.l2| ). This gives 



F" - 2A'F' - AA"F 



2^F' + 4A'^F 

00 00 



nF, 



(3.17) 



to be solved in the bulk. This is the principle equation that will be studied and solved 
below. We note in passing that each eigenmode aF„ = — m^F„ to this equation has two 
integration constants and one mass eigenvalue. One constant corresponds to the overall 
normalization. The remaining integration constant is fixed by the boundary condition at the 
Planck brane, and the mass is determined by the boundary condition on the TeV brane. In 
the stiff potential approximation we use the boundary condition given by ( p.l5| ). 

It is then possible to show that a solution F to the above equation automatically implies 
that the equation and the remaining Einstein equation are satisfied. In particular, starting 
with ( p.l6|) , one uses the derivative of (|3.12| ) to eliminate F". The resulting equation, call 
it E, is then differentiated and the combination = E' — 2A'E is constructed. Using the 
background field equations and (|3.12| ) one arrives at the ip equation 
equation is obtained from the equation after substituting for ip'. 



Finally, the 6R 



4 General Properties of the Equation 



First we show that the single ordinary differential equation for F(y) given in ( |3.17| ) can 
always be brought into the Schrodinger form. For this we first transform the equations into 
the coordinate system where the background metric is conformally flat. This is achieved 
by the change of variables dze~^^^^ = dy, where A{z) = A{y{z)). In these coordinates the 
equation simplifies to 

F" - 3A'F' - AA'F - 2^-F' + A^-AF = -m^F. (4.1) 

00 00 

After the rescaling of the field F hj F = e^^'^^cp^F we obtain the Schrodinger-like equation 
-F"+{ ^A'^ + ^A" -A'§ + 2(§X -^\F = m'F. (4.2) 



However, this by itself does not guarantee hermiticity of the differential operator in ( |4.2|) . 
The reason is that this operator is defined only on a finite strip, and therefore in addition 
to writing the equation in a Schrodinger form one also has to ensure that one has hermitian 



8 



boundary conditions for F. For the differential operator in ( ^.2| ) to actually be hermitian on 
the finite strip between the two branes, one also has to require that for any two functions 
Fi, F2 on the strip F{(0)F2(0) -Fi(0)F^(0) -Fi'(zfe)F2(2b) = 0, where and 

denote the positions of the branes in the conformally flat z coordinates. Once this condition 
is satisfied, it is automatically guaranteed by the usual theorems that all eigenvalues are 
real, that the eigenf unctions are orthogonal to each other and that they form a complete 
set. The actual boundary conditions that F has to satisfy can be derived from the general 
boundary condition given in ( |3.14 ). In the particular model considered in this paper the 



boundary condition in the y coordinates is given by 

± = 72 ^ ± 2u(t)±F (4.3) 

In the special limit when 7-1- ^ 00 this boundary condition reduces to <y9 = on the two 
boundaries, which together with the constraint equation ( p.l2| ) between if and F just implies 

(F'-2v4'F)|i = (4.4) 

at the two branes. Upon transforming to the Schrodinger basis and z coordinates the bound- 
ary condition will be replaced by 

F' = F{]-A' - f ) (4.5) 
z (Pq 

at the branes. This boundary condition clearly satisfies the hermiticity properties and thus 
will ensure the appearance of only real mass eigenvalues of the coupled system. This will 
also be the case that we will analyze in full detail in the following sections. As for the general 
case, when 7^ is finite, the boundary condition will not be hermitian. This can be easily seen 
from the fact that the general boundary condition involves y?' at the branes, which should be 
expressed from (|3.12|) in terms of F", F' and F at the brane. The appearance of F" in the 



boundary condition will generically ruin the hermiticity of the operator. Nevertheless, one 
may eliminate F in favor of the eigenvalue, and one can in principle solve for F. The non- 
hermiticity by itself however does not mean that the eigenvalues are not real. In fact, since 
is a real scalar, and F a component of the metric tensor, both of these functions have to be 
real to start with, which guarantees at least the appearance of only real eigenvalues. While 
for the model studied here (see Section ^ the radion is not tachyonic, it is unclear whether 
for a general potential this remains true. However, the orthogonality of the solutions is not 
guaranteed by anything, and will likely be violated in general for the non-hermitian boundary 
conditions. It would be interesting to understand the physics behind the non-orthogonality 
of these solutions in more detail. 



5 Approximate solution for the KK tower 

We have seen that the coupled radion-scalar system leads to a single ordinary second order 
differential equation. From now on we will always assume that we can use the limit 7^ — 00, 
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and be able to use the hermitian boundary conditions ( [4.41 ). In the following we will present 
an approximate solution to these equations. For this, we will first neglect the backreaction 
of the non-vanishing scalar background on the metric. This will lead us to a simple Bessel- 
type equation, which will give a very good approximation for the masses of the KK tower 
of the fields. However, surprisingly, in this approximation the radion (which we identify 
as the lowest lying solution of (|3.17| ) remains massless. Therefore, after presenting this 
approximation, we will give a perturbative analysis for the effect of the backreaction of the 
metric on the radion mass. We will find that as expected, the radion mass will be of order 
TeV, but somewhat lighter just as predicted in [15[| . 

To find the actual wave functions and masses for the radion-scalar system, we will use 
the particular model put forward by de Wolfe et al. ||2^ and summarized in Section ^. First 



we neglect the backreaction of the scalar field background on the metric, which seems to be 
a good approximation as long as K,(f)p^T -C 1. In this case the equation for the radion field F 
reduces in the Schrodinger frame to: 

where a is given by 

„ = ----. (5.2) 

In these coordinates the boundary conditions at the brane simplify to 

„, Oik _ 

F' + -F = 5.3 

kz + 1 ^ ' 

at the locations of the branes at z = and Zb = -^(e^^'" — 1). The solutions of these equations 
are given by linear combinations of the Bessel functions Ja+i and Neumann functions N^:i: 



1 1 1 1 1 

Fn{z) = an{z + -)^N^^{m^{z + -)) + 6„(z + J^^{m^iz + ^)). (5.4) 



The mass eigenvalues m„ can be determined from the boundary condition ( |5.3| ). Using the 
relation for Bessel functions 

Z'^{x) = Zn-^l{x) Zn{x) (5.5) 

X 

the boundary conditions at the two branes simply reduce to 

aniV„_.(^) + 6„4_.(^) = 0, 

a„Ar„_i( — - — ) + fe„4_i( — ^ — )=0, (5.6) 
which yields the simple equation 

K-)^^„-,(^)^^i5S^-4-i(^)^o, (5.r) 
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Figure 1: The lowest mass eigenvalues for the coupled radion-scalar system for u/k = 1 
are given by the zeroes of the function b{m) defined in (|5.7| ). On this plot m is given in 
units ke~^'^°, therefore the mass spacings are given by the TeV scale. Note, that in the 
approximation leading to this equation the lowest lying state is still massless. 



which can be used to determine the mass eigenvalues This can be done numerically. In 
Fig. |I|we show the lowest mass eigenvalues for a = —2.5, which corresponds to the somewhat 
unrealistic value | = 1. In Fig. ^ we show the dependence of the first non- vanishing mass 
eigenvalue on the value of a = —3/2 — u/k. One can easily see from ( |5.7| ), that m = is 
always a solution to ( |5.7| ), therefore in the approximation we are using the radion is still 
massless. For the higher states of the KK tower it is a good approximation to use the mass 
eigenvalues obtained from ( ^.7| ), because the masses are of the order (and even larger) than 
the TeV scale, thus in the limit of small backreaction that we are considering throughout the 
paper these masses will be only slightly modified. The radion (which appeared as the zero 
mode above) however needs special treatment, since the shift in the mass (which is usually 
negligible for the higher KK modes) coming from the backreaction of the metric background 
due to the scalar field is the leading order contribution to the mass for the radion. Below we 
will estimate the size of the radion mass in perturbation theory. 



6 The radion mass 

In the previous section we have seen what the approximate wave-functions and masses are 
for the KK tower of the coupled radion-scalar system. In this approximation of neglecting 
the backreaction, however, we have still found a vanishing radion mass. This is in fact easy 
to show for a general stabilizing potential. From (|3.17| ), F = e^^ is always a solution with 



zero mass if A" is neglected in the bulk. Thus the radion mass is always proportional to 
the backreaction of the metric independently of the details of the potential of the stabilizing 
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Figure 2: The dependence of the mass of the first KK mode on a. Here mi is again given in 
units ke~^'^° and is therefore of the order of the TeV scale. 



scalar field. In the following, we will show how the backreaction generates a non-vanishing 
mass for the radion field. For this, we start with the equation describing the radion wave 
function in the y coordinates: 

F" - 2AF' - AA'F + 2uF' - AuA'F + m^e^'^F = 0, (6.1) 

where A{y) is given in ( |2.19| ). The appropriate boundary condition is F' — 2A'F = at the 
branes. In the special limit 7-1- ^ oo the other boundary condition is = 0. Thus we will 
treat the backreaction as a perturbation, and look for the solution in terms of a perturbative 
series in I = k0p/v^. Then we write the solution as 

Fo = e2^l^l(l + I'Uy)) , ml = I'm' , A{y) = k\y\ + ^e-^^I^L (6.2) 

6 

Expanding the solution as above and keeping only the leading terms in P we obtain the 
equation 

+ 2{k + u)fl = -m2g2fc|j;| _ _ u)ue-'''\y\ (6.3) 
along with the boundary conditions 

/o + \ue-'^\y\ = (6.4) 

at the location of the branes. One can easily find the most general solution for / from the 
equation in the bulk, which is given by 

m = Ce-'^^^-^\y\ - ^(^^e^^l^' - ^^^^e--l^l, (6.5) 

12 



where the integration constant C along with the radion mass m is determined by the bound- 
ary conditions at the brane. This way we obtain the radion mass to be 

''''radion ~ ' ^ 

where tq denotes the location of the brane. Note that this result is very similar to the answer 
obtained from the effective theory computation using the naive ansatzM, IIH], except for the 



important difference in the power of u/k. The exact result obtained here scales as (u/k)^, 
whereas the effective theory result would scale as (M/fc)3/2.[] It would be very interesting to 
understand the origin of this different scaling. For this model to give the correct value of the 
weak scale without reintroducing a large fine-tuning again one needs u/k k, 1/37, thus the 
radion mass turns out to be somewhat lighter than the TeV scale. It is suppressed by the 
factor /|e~"''° compared to the TeV scale. Thus in this approximation rriradion ~ ^ TeV, 
which could be at least in the range of a few GeV's. Of course, we need to emphasize that / 
is not necessarily small for the stabilization mechanism to work, we took this limit only for 
calculational convenience. 



7 Coupling to SM fields 

In this section the coupling of the radion and KK tower of to the TeV brane are obtained. 
In particular we demonstrate that the bulk scalar field gives a small correction to the radion 
kinetic term, and thus the kinetic terms obtained from the Einstein-Hilbert part of the 



action dominate, justifying the results obtained using the naive ansatz |1J, |T5[. 

In the previous section it was seen that by including the backreaction an (9(TeV^) for 
the radion is obtained. The wavefunction is then 



Fo(x,y)=e2%l(l + /Vo(z/))i?(x) (7.1) 



where foiy) is given by the integral of ( p.5|) . Since the radion mass is O(TeV^), and by 
assumption -C 1, we see by inspection that the backreaction induces a small correction 
to the unperturbed wavefunction. So for the purposes of determining the coupling of the 
radion to the TeV brane it is sufficient to include only the unperturbed wavefunction, namely 
F(x,y) = e^'^'^'i?(a:). Then a straightforward calculation gives 

-M^ J dy^n D GM^dR)'' I e-^^e^'^l^l = ^{e^'^^o - l){dRy (7.2) 

So the normalized radion r{x) is R{x) = r{x)e-'''''>/V6Mpi, since M^/k = Mpi/2. This 
implies a coupling to the TeV brane fields which is 

^2fcro / 1 , n(l'^\ \ Tt^T^ — ^^^^ TvT^ (^ _L r)('/2- 



i?(x)e^'=^« (1 + 0{l')) TrT,. = y^^^^^^TrT,, (l + 0{l')) , (7.3) 



*We thank Jim Clinc for these observations. 
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where the left hand side of this equation is a consequence of the fact that the induced metric 
on the TeV brane is e~^^''^°''(l — e'^^^°R{x)'q^y). The couphng obtained this way agrees 
precisely with ( p.4|) . This is perhaps surprising, since the latter computation used an ansatz 
which did not satisfy the equations of motion. This makes us suspect that results which 
depend only on the leading unperturbed form of the radion wave function will be correctly 
captured by the naive ansatz. 

Now we address the issue that was originally raised by [^. Is the radion kinetic term 
dominated by the kinetic term of the bulk scalar field or the bulk gravity action, and in 
particular is the former hierarchically larger? To answer this, we need the change in (p 
caused by a fluctuation in the radion. Since ip = when the backreaction is not included, we 
must include the leading backreaction correction to the radion wavefunction, given by the 
integral of ( |6.5| ). From ( p.l2| ) we compute that the change in ip to 0{P) due to the radion is 



where 



My) ^ /^(y) + ^ue-2«M 

2{2k + u) 3k ^ ' 

This fluctuation in (f) then contributes to the radion kinetic term at 0(/^) an amount 

/q/2 n 
dye-'^g^^d,vd,^ = ^{dR)' j dye^^''-^\y\ h{yf ■ (7-6) 

From ( [7.2|) the unnormalized contribution from bulk gravity to the kinetic term is ~ e^'^'~°. 
So we only need to consider those contributions from which are comparable or larger to 
this. Recalling that m? ~ e~'^^'^° , it is seen that the largest terms in ( [7.6| ) are at best e^'^''°. 
Explicitly performing the integral one finds that it is 

5C = 2l'4r.m?e''^^-'-- (-^ + - (l + 0{f)) . (7.7) 

K^/e^ \6k — u k — 6U k — uj^ ' 



This is typically ~ t^v?(?^'^^M'^ jk^ ^ which is smaller than (|7.2| ) since we assuming that the 



backreaction is small, I <^ 1, and also that u <^ A; to obtain a realistic hierarchy. So the 
radion kinetic term is dominated by the contribution from the bulk gravity, and receives a 
small correction from the stabilizing bulk scalar field. 

In Section ^ it was found that for the simple boundary conditions = (corresponding 
to the limit d'^\±/d(j)'^ 3> 1) a self-adjoint equation for F was obtained. The general solution 
to this is 

F{x,y) = Y.anFn{x,y) (7.8) 
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where F„ is a mass eigenstate, and the a^s are some numbers. We expect that F includes the 
massive radion, but where did all the other states come from? It is helpful to reconsider what 
happens when the backreaction is neglected. In this limit the KK tower in F completely 
disappears and only the (massless) radion remains. This may be observed from (|3.12| ), 
since neglecting the backreaction corresponds to k^^q ^ 1, and this amounts to setting 
F' — 2A'F = 0. The only solution for F in this case is the radion zero mode F = e^'^'^L 
Once the backreaction is included, however, the fluctuating modes in and F are correlated 
through ( p. 12 ). In particular, a general fluctuation ip induces a change in F. The sum over 



KK states appearing above is then just the decomposition of F into these KK eigenstates. 
It is then expected that the coefficients for the non-radion states to be suppressed by the 
backreaction. 

The preceding remarks imply that the TeV brane fields, which couple to the induced 
metric F, also directly couple to the KK tower, by an amount suppressed by the backre- 
action.Q Since F ~ (p'^ip is already suppressed by the backreaction, therefore in order to 
compute the induced metric to lowest order in the backreaction we can use the zeroth order 
wavefunctions for ip. The normalized KK fields are given by 

ifnix, Z) = ^^{kz + lY-h+u/kim^Z + l/k) (l + /2/n(l/)) • (7.9) 

Here ipnix) are the normalized 4-D fields satisfying Oipn = —T^n'^n- The orthogonality of 
these solutions when the backreaction vanishes (/ = 0) follows from the boundary condition 
(fin = and the properties of the Bessel functions. Also ipn{x) are the normalized 4-D fields 
satisfying = —m^ipn- The normalization constant is 



iV„ = 4=e'='-oj3+„/4!^^e'='^»). (7.10) 



As discussed previously, the lowest order masses m„ are determined by J2+u/fc(e^'^''°m„//c) = 
and are real since the operator equation with these boundary conditions is self-adjoint. The 
coupling of these fields to the TeV brane is given by 

F(x,?/ = ro)TrT^, (7.11) 

where F is the solution of (|3.12| ) for the solutions given above. One finds the coupling 

^ -A,x>„(a;)TrT^,. (7.12) 



rrir. 



The model-dependent couplings that appear are 



\s = ^uK^e''^''' ~ O(^) (7.13) 
3 k 



*The KK modes of the scalar field do not mix with the KK modes of the graviton. The reason is that the 
only way they could mix is by a coupling of the 5D trace of the metric to the scalar KK modes. However, the 
graviton is traceless, and the trace is basically identified with the radion, therefore no additional graviton- 
scalar mixing could be introduced. 
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and 

is a numerical constant of 0(1). While the inclusion of the backreaction leads to a TeV 
suppressed coupling for the KK modes, the size decreases rather rapidly due to the l/m„ ~ 
1/TeV suppression, as may be observed from inspecting Fig. 1. 

The coupling discussed here implies that the KK modes of (p can be directly produced 
at future colliders, and they also decay directly to Standard Model fields. This may be 
puzzling at first, since the stabilizing potential may have a global discrete symmetry, such 
as Z2, which would naively imply that some of these KK modes are stable. The background 
vev for explicitly breaks this symmetry, however, and this allows for all the KK modes to 
decay into the brane world fields. 

The direct coupling of the KK modes from the stabilizing fields may have interesting 
implications for search strategies and current limits on the Randall- Sundrum framework. In 
particular, it may be important to not neglect the stabilizing potential when discussing these 
issues. However, when the backreaction is small, the size of their couplings is suppressed by 
u/k ~ 1/40 compared to the that of the radion. Therefore, in what follows, we neglect these 
states in the loop computations. 



8 Cosmological Implications 



The subject of brane cosmology has recently attracted lots of interest [|10|, |Il|, |1J, |31|, ^ 



p3| , [3^ , p5| , |36| , p7| . Most of this was due to the realization, that the expansion of a brane 
universe could be significantly different from the ordinary Friedmann-Robertson-Walker 
(FRW) cosmology |T^, |TT|. However, it did not take very long to realize, that this is simply 
due to the fact, that a generic brane model (like the one presented in |]10|) can not give the 
ordinary cosmological evolution, since gravity is in general manifestly higher dimensional. 
This means that in these models the 4D effective theory is usually not described by ordinary 
Einstein gravity, but generically a complicated scalar-tensor theory of gravity. However, 
observations show that our Universe is described by Einstein's theory of relativity to a high 
precision, therefore one has to require from the outset that a brane model reproduces ordinary 
Einstein gravity, at least at long enough distances. Once this is achieved, the cosmological 
expansion will be automatically described by the ordinary FRW Universe, which simply 
follows from the fact that the effective theory is ordinary Einstein gravity. Thus one can 
see that the issue of unconventional cosmologies is nothing else but the issue of whether one 
recovers 4D gravity. This issue manifests itself in the case of the Randall-Sundrum two- 
brane model due to the presence of the radion field. Without a stabilizing potential, the 
radion field will be massless, and yield additional long range forces, and also contribute to 
the expansion of the Universe, yielding an unconventional cosmology, which is presumably 



excluded by the requirement for a successful nucleosynthesis [|TT|. Thus the radion field has 
to obtain a mass. Once it is massive, gravity on both branes will be ordinary 4D gravity, 
and thus the cosmology will be conventional below temperatures comparable to the radion 
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mass. This has been explained in great detail in |T^, and also in |31|. In [|I^ a simplified 
calculation for the cosmological expansion has been presented, where the wave function of the 
radion has been neglected, and also the effects of the stabilizing scalar field were included 
by adding a five dimensional potential for the radion field V{b). Assuming that that the 
potential V{b) is very steep, it was shown from a perturbative solution of the bulk equations 
that the ordinary FRW Universe is recovered. It was also argued, that the 55 component 
of the Einstein equation, which in the absence of a stabilizing potential usually leads to the 
unconventional expansion equations will only determine the shift in the radion field due to 
matter on the wall, and does not result in unconventional cosmologies once the radius is 



stabilized. Below we demonstrate, that the results of [|I4| which were neglecting the radion 
wave function, and also did not include the fiuctuations of the scalar field at the brane 
remain valid in the more precise framework of radion stabilization explained in the previous 
sections. In particular, we will show, that the result obtained for the shift in the radion field 
due to matter on the walls in Eq. (4.15) of [11^ | is exactly reproduced in the full calculation. 
To compute G55 we use the ansatz 



ds'^ = n(t, uYdt^ — a^{t, y)d^x — h{t, y^dy^ 



for which 



G 



55 




a / a 
a \a 



(8.2) 



The jump equations for a and n on the TeV brane imply 



a 



-y(A-(0) + p)& 



n 



K 



:-A_(0) + 3p + 2p)6 



(8.3) 



Here p and p are the bare energy matter density and pressure on the TeV brane, which are 
related to the physically measured quantities on the TeV brane by po = pe"'^'^'', etc. where 
e~^° is the scale factor on the TeV brane. Then averaging the G55 equation about the TeV 
brane and linearizing to 0(p, -F, ^p) gives 



< Gss > 



,A2 




-(3P-P) 



^.4) 



Terms with h are higher order in p and are dropped, and 6 = 1 + 2F has been used. For 
late-time cosmology in the presence of radion stabilization it is reasonable to use the FRW 
equation 



+ 



Oo 

ao 



1 



(3p^ - p* + 3po - Po) • 



^.5) 



Note this includes a contribution from matter (p^,) on the Planck brane. Also implicit in 
the use of this equation is the assumption that the time variation of the radion is negligible. 
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which is justified a posteriori. Then using the relation k^A_(0o) = —6(1 — e '^"^°)/Mpi gives 

<G55> = ^^^4^ + ^(3Po-Po + e-^^"(3p*-p*)) 

^ /.^A-(0o)aA_(0o) ^^ ^ 2^^A^_(0o) ^^ ^g g^ 



The G55 equation is 

G55 = f^'n, = (^0'^ - ^755(^(V0)2 - V(0))) . (8.7) 
Then the averaging of T55 and hnearizing using (^.2]) gives 



/t^ < T55 >= - V{(Po)) + ^\4>'oV' - ^FV - —^) (8.8) 

with all quantities are evaluated on the TeV brane. Using the background bulk equation 
( p.lOp and the jump equation ( |2.12| ) the leading terms are seen to cancel. Then after using 
the background equations, the jump equations for the background fields, and some algebra 
gives 

g4^o , \ 2 dV 

^ (3po - Po + e-2^o(3p. - p.)) = k2(0'o<^' _ 20 - — <^ - 4A>[,<^). (8.9) 

Using ( p. 16 ) to eliminate y^', ( 3.17 ) to eliminate F" in favor of the mass eigenvalue, and 
( p.ll| ) to eliminate 0q finally gives 

^ (3po - Po + e-2^"(3p, - p,)) = -3e2^0m2F. (8.10) 

But the shift in the distance between the two branes is obtained by integrating the line 
element, which gives dr^ = i?(e^^^'' — l)/k. Then since F = Re^^'^^\ one obtains 

which is precisely the result found in [|14| obtained by using a 4D effective theory.^ This is 
perhaps not surprising, since for constant radion field the naive ansatz and the full metric 
including the wavefunction of the radion are equivalent up to a coordinate transformation. 
So in an adiabatic approximation the leading order result using the naive ansatz should 
agree with that obtained from using the correct radion wavefunction, if the fiuctuations in 
the scalar field are ignored. It is less clear why the full answer including the contribution 
from the scalar field turns out to be exactly equal to the calculation using the naive ansatz. 
Note that matter on the Planck brane causes a smaller shift in the radion compared to an 
equal amount of matter on the TeV brane. This is because the radion wavefunction is peaked 
at the TeV brane, and it couples more weakly to the Planck brane relative to the TeV brane 
by precisely the amount e"^"^". Thus one finds the very general result that in the presence 
of matter on the branes and a stabilizing mechanism, the 6*55 equation determines the shift 
in the radion. 

*A translation dictionary between two different notations is required: hr^ — mQbQ/2. 
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9 The Effective 4D Lagrangian 



In the previous sections we have argued that in the presence of a stabihzing potential the 
hnear couphngs of the radion and bulk scalars is given by 



\{drf - i^^^nf - ml^t) + DH^DH 

+ + E TrT,. + iH^Hn - V{H). (9.1) 

The masses appearing here are O(TeV), and their particular value depends on the details of 
the stabilizing mechanism. The scale A = e~^'^°Mpi in the Randall- Sundrum model, but here 
we have left it general. The other scales are A„ ~ m„, and the a„ are also model-dependent, 
and vanish in the limit of small backreaction. In the remaining sections we restrict ourselves 
to the above Lagrangian, and do not commit ourselves to any specific mechanism of radius 
stabilization. For the electroweak analysis we neglect the contributions of the KK modes 
from (f). 

Note that in the above Lagrangian we have also included a curvature Higgs scalar operator 
WHTZ. The presence of this operator leads to interesting signals for discovering the Higgs 
and radion at future colliders [^. In particular, the branching fractions of the Higgs and 
radion to gg and bb can be substantial different from that of the SM Higgs. 



As discussed in the presence of the conformal term WHR leads to both kinetic 
and mass mixing between the neutral Higgs and radion. Below we summarize the relevant 
formulae for mixing and couplings. The interested reader is referred to the next section for 
details. 

One finds that the "gauge" h and r are related to the mass eigenstates hm and by 

h = (cos 9 — sin 9)hm + (sin 9 H — — cos 9)rm , (9.2) 



where 



r = cos6'^ — sin6'^ , (9.3) 

Z Zj 

2 

tan 29 = U^jZ— , , (9.4) 

m;. — mf^ — 647^(1 — I24) 

7 = -^,Z2 = 1 + 6^7(1 -60, (9.5) 

where v ~ 246 GeV is the electroweak vev. Requiring that the quantity be positive (in 
order to avoid ghost-hke states) places an upper bound on the value of ^, for a given 7. 
Physically this requirement comes from maintaining positive definite kinetic terms for h and 
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In this basis, the couphngs of the physical radion and Higgs appropriate for tree-level 
studies are 

- ({cose- (6e - l)^^)hm + (sine + (6e - l)^^)r„^ TrT^.- (9.6) 
In the ^ ^ limit one recovers 

-{h- 7r) TrT^, , (9.7) 
obtained in [jn|, |15|. Note that TrT^,^ includes SM Higgs contributions. 



10 Curvature— Scalar Mixing 

In this Section the effect of introducing a curvature-scalar interaction are reviewed. The 
discussion parallels |^|, however some of the resulting formulae are slightly different because 



here terms of 0(7^) and 0(7^^^) are kept. 

We begin with the couplings of the radion and Higgs to the SM fields before electroweak 
symmetry breaking. The induced metric on the TeV wall is 

^-/(x) = e-^^(-)-^^^^-«(^')^,.(x), (10.1) 

where the warp factor includes the backreaction, although its inclusion is not necessary for 
our purposes. The canonically normalized radion r is 

So we express the induced metric expanded about a Minkowski metric as 

= e-^^^^'^-^'-<^\, ^ e~^^^^''^^\r)r^,, (10.3) 

with 

7 = ^,A = M„e-'='^». (10.4) 

The four dimensional effective action we consider is 

STeV = j d^x^a(gZdD,H^D,H -V{H)) 

+ j d^x^^ ((Vr)2 - my) + J d'x^,^R{g,,,,)H^H + Ssm- (10.5) 

To canonically normalize the Higgs and other SM fields, we perform the field-independent 
redefinition 

H ^ e'^^'^^H , ijj e=^^(''«)/V (10-6) 
In this basis the Higgs-radion potential is 

V{H,r) = n\r)V{H). (10.7) 
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Note that V also includes the effective 4-dimensional cosmological constant, which we assume 
to vanish. Clearly this potential has a minimum at the same location as V{H), so that the 
EWSB vacuum is r = and = v/V2. 

We consider the presence of the curvature mixing term 

A = V9^dW9ind)H^H. (10.8) 

Our choice of signs for ^ is such that the Higgs potential receives a positive mass-squared 
correction in a de Sitter phase when ^ is positive. Since this is a renormalizable interaction, 
there is no reason for it not to be present, or to be suppressed. What makes this operator 
important in this case is that TZ contains the induced metric, rather than just the ordinary 
4-dimensional metric. In particular, 

n(n^(r)r]^^) = -6Q-^ (nln Q + (V In Q)^) . (10.9) 

So the curvature-scalar interaction is 

jC^ = -6^Q2 (^QinQ + (VlnQ)^) H^fH (10.10) 

To see the effect of the curvature scalar interaction we expand = {v + h)/ \pl and i7(r) = 
1 — 7r/t' + • • •. We need to only expand f2 to linear order since the derivative terms are 
already of 0{r). This gives at quadratic order 

= + "ii-i^idrf. (10.11) 

where a total derivative has been dropped. The ^ terms clearly introduce kinetic mixing. 
The full radion-Higgs Lagrangian to be diagonahzed is 

111 1 
£ = — /o/i - -m\h^ - -(1 + 6^72)rnr - -m^r^ + (10.12) 
2 2 2 2 

The mass parameters m^, ruh are the masses of the radion and Higgs, respectively, in the 
limit ^ = 0. The kinetic terms are diagonahzed by the shift h — h' -\-&^'yr' / Z , and r = r' /Z. 
Here 

= 1 + 6^7^(1 -60 (10.13) 

is the coefficient of the radion kinetic term after undoing the kinetic mixing, and is therefore 
required to be positive in order to keep the radion kinetic term positive definite. For a fixed 
cutoff A this restricts the size of the mixing parameter ^. It must lie in the range 




(10.14) 



for non-zero values of 7. Otherwise one has a ghost-like radion field, which presumably 
signals an instability of the theory. 
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This rescaling diagonalizes the kinetic terms, but introduces mixing in the mass matrix. 
A final rotation h' = cos 6hm + sin and r' = cos 9rm — sin 9h.,n brings the Lagrangian to 
canonical form. With the above definition of the sign of the rotation, the rotation angle is 

2 

tan2^ = 12^7^^ J^^ (10.15) 

We note that for moderate values of ^ and 7 (i.e., Z"^ > 36^^7^) the mixing angle tan 2^ is 
negative when > m^. For small 7 we can expand 

2 

tan 2^ - 12^7 J^^ , + 0(7"). (10.16) 
m; — 

Putting everything together, the relation between the gauge and mass eigenstates is 

h = (cos^-^sin^)/i^ + (sin^ + ^cose)r^ , (10.17) 

r = cos^^-sin^%. (10.18) 



The mass eigenvalues are easily obtained 
1 

2^2 



m| = (ml + (1 + 6^7')m^ ± (K - + 6^7^))' + \^^l\^mlmlfl'') (10.19) 



The heavier state (+) is identified with the state with the larger of {ml^my). 



11 Radion Couplings and Feynman Rules 

In this section we derive the Feynman rules relevant to the computation of the oblique 

parameters S. T, U . 

Before proceeding, we pause to ask whether higher-order couplings such as 

(11.1) 

also affect in particular the electroweak precision measurements S, T and U. This operator 
could either be directly present, or generated from the above linear coupling due to a non- 
trivial kinetic term for the radion. Although this operator contributes at one-loop to the 
gauge boson two point functions, it is easy to sec that they do not contribute to T, since the 
my contained in TrT^j, is canceled by the 1/my appearing in the expression for T, nor to 
S or U since the contribution of this operator to the vacuum polarizations is momentum- 
independent. Thus we need to only consider the linear coupling 

—TrT. (11.2) 

V 

This operator will have a contribution to the oblique corrections similar to that of the 
standard model Higgs. First we discuss the Feynman rules for the interactions from the 
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^TrT operator. The interaction Lagrangian term relevant for the gauge-boson propagator 
corrections is just given by 

= - Jr(2M^W;iy^- + MlZ^Z^^). (11.3) 

In addition, to ensure gauge invariance of the results, one also has to examine the gauge 
fixing terms carefully. The gauge fixing Lagrangian for the W and the Z are given by 



a la 



(11- , 

where the \l/'s are the would-be-Goldstone bosons, and a is the gauge fixing parameter 
in the gauge. Note, that since the gravitational background is non-trivial, we have 
replaced the ordinary derivatives by covariant derivatives. The background metric is given 
by g^iv = ^'^{f^)Vfj.u = ^~~^fivi therefore the covariant derivative of a vector will take the 
form 

D^V^ = n-\d^V^ - '^d^rV^). (11.5) 

Thus, from the gauge fixing terms one also obtains three-point interaction vertices of the 
form 

-^9/ = —{d^Z>'){d,rZ') + ^(a^M/A^+)(9,rW^^-) + ^(a^M^^-)(5,rW^^+). (11.6) 
va va va 

With these operators added, the Feynman rule for the three-point function is given by 

Vtiu {P2f,Plu + P3uPlf,) (11-7) 

In addition to the cubic vertices evaluated above, there are also 4-point couplings of the 
radion and the gauge bosons. These terms will not contribute to the oblique electroweak 
corrections, however they will be important to obtain a gauge invariant answer for the gauge 
boson vacuum polarization diagrams. These terms arise from two different sources. The first 
source is the conformal coupling — e~~TrT to the trace of the energy momentum tensor. 
In the formalism of |T^ this can be obtained by a very careful expansion of the interaction 
terms + -^^j^dr[W DH + h.c). Either way one finds the additional operators 

±.r\^W;W^- + MlZ^Z^^). (11.8) 

The other source of quartic interaction terms are again the gauge fixing terms. One simply 
expands these to higher order to obtain the interaction terms 

27^ 



d^rd^riZ^Z" + 2W''+W''-). (11.9) 
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The interaction terms involving the would-be Goldstone bosons just contribute a total deriva- 
tive, and thus they can be omitted. The above operators give rise to the following Feynman 
rule for the 4 point function: 



5 Vl^u H ^[Pl^iP2v + PlvP2^i) (11.10) 




The presence of the terms proportional to ^ in ( p. 1 . 7| ) and ( p.l.lO| ) are in fact crucially 
important to obtain gauge invariant amplitudes, however for calculations in the unitary 
gauge a ^ oo their effect vanishes. 

This is however not the complete story for the Feynman rules. The reason is that the 
radion couples conformally to the metric and so is not the same as the Higgs. Thus for 
example one finds that at one-loop the radion has the anomalous coupling 

X&G^G^.G'^^ (11.11) 

in addition to the usual momentum-dependent coupling obtained from one-loop diagrams 
with internal fermions. Here he is the beta-function. This may be understood as due to the 
scaling anomaly together with r as a generator of scale transformations. Diagrammatically 
this result is obtained by preserving the conformal coupling of r when the theory is regulated. 
For dimensional regularization this means that the radion must couple conformally to the 
D-dimensional metric. Since the linear coupling of the radion is obtained from varying the 
induced metric, for loop computations the radion should couple instead to Tr^T^,^, but where 
now the trace is evaluated in D-dimensions. This differs from the above coupling by some 
operators whose coefficient vanishes when D ^ A. Since this e = A — D dependence can 
be offset by poles appearing in the loops, the appearance of these additional operators can 
result in finite non-zero results in the D — > 4 limit. These operators will indeed have a 
non-vanishing contribution to the S and U parameter. Next we calculate the Feynman rules 
for these "anomalous" couplings. The interactions we should therefore study are 

C = --TiiT + —TtdT. (11.12) 

V V 

We point out that it is the original 'gauge' radion, r, that has the conformal coupling to 
the metric, and consequently it is this field which appears in the above interactions. But 
a result of the curvature-scalar term interaction is to introduce mixing between the radion 
and Higgs, which implies that after transforming to the mass basis the physical Higgs, hm, 
will have couplings similar to those above. 

The radion will thus have the interaction terms {D = A — e) 

L = j-TtdT^.u = - -f^-F,,F^'' + ^^-{2M^W+W- + MlZ'), (11.13) 
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where the first term is the one we have already discussed above. We will later add in the 
Higgs and mixing coefficients. The terms relevant to the gauge boson propagators from the 
last two terms are: 



Qanom 



€7 

2M^W^W''- - MlZ^Z" 



;il.l4) 



In D dimensions, one also has to modify the gauge fixing terms. The covariant derivative of 
a vector field will be modified to 



;ii.i5) 



In addition, the y/g factor in front of the gauge fixing terms have to be modified to . 
Thus these interaction terms modify the Feynman rules for the interaction vertex given in 
( ITTl ) to 



( 1 - - ) r]^^ H (P2 ■ PsVf^u - P2uP3f, ) {P2fMPlu +P3uPlf,) P2fMP3u ■ 11.16) 

V 2 V av av 

In addition, the four-point vertex in (|11.10|) is also modified by terms proportional to e, 
which however do not contribute to a calculation in the unitary gauge. 

We close this section by discussing how to take the mixing between the Higgs and the 
radion due to the possible presence of the curvature-scalar mixing operator into account. 
The interaction in the gauge basis is 



Using 



C = Tr4T + — Tr^T. 

V V 



h = ahm + br^ , 



;ii.i7) 



;ii.i8) 

11.19) 



where the coefficients a, b, c and d can be read off from ( |10.17| ),( |10.18| ), then 



-(a - 7c)— + hd - b) — 

V V 



Tr.T 



11.20) 



Thus the Feynman rules for the radion (and also for the Higgs) have to be modified such, that 
the above mixing terms are taken properly into account, for example for the mass eigenstate 
radion the Feynman rule will be 



i^de i{2-e)^d i^ed 

m,^v + — (P2 -PS^^M^ -P2uP3p,) — [P2t,Plv +P3uPlf,) ' 

11.21) 
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12 Electroweak Precision Measurements 



In this Section we consider the corrections of the Randall-Sundrum model to the obhque 
parameters. Our analysis also applies more generally to a model with a light scalar coupled 
conformally to the SM metric but with a typical coupling of 0{TeV)~^. 

Our approach is to use an effective theory with cutoff 0{A), similar to the approach taken 
in |25|. Below this scale the only light fields are the radion and Higgs whose contributions 
we are going to calculate explicitly. In our approach the effect of any modes heavier than the 
cutoff are included by introducing higher dimension operators that directly contribute to the 
oblique parameters. This in principle includes the effects of the heavy spin-2 KK states, for 
example, which are typically heavier than the radion. A direct computation of the effect of 
the heavy spin-2 states using a momentum-dependent regulator has been presented in 



In the previous Section the radion coupling to the gauge bosons was obtained and found 
to be similar to that of the Higgs. The contribution of the Higgs to the oblique parameters is 
by itself divergent, but these divergences are canceled by the contribution from the pseudo- 
Goldstone bosons (or the longitudinal states of the massive gauge bosons). Thus for the 
radion one expects a divergent contribution, but in contrast to the Higgs there is no additional 
source to cancel this. This is perhaps not surprising since the radion interactions are non- 
renormalizable. 

A set of operators that provide the necessary counterterms for the wavefunction renor- 
malization is 

(h^HTtW^.W'"' + I- tan^ OwH^ H B^^.B^^" + tan OwH^W^.B^" h] , (12. 



Ox = '-jf- [H^HTiW^.W^'' + - tan^ OwH^ H B^.^B^'^ + tanOwH^W^^B'^'' H ^ , (12.1) 

where Wf^u = W"'t"', with the generators normalized to 1/2. Note that the last operator is 
gauge invariant, since for a gauge transformation U, W^u — ^ UW^i^W . Setting the Higgs to 
its vev in the above operator gives 

fx^^(w'i«'-"" + ^^.^^-). (1") 
In this model the radion does not contribute to 77 and 7Z wavefunction renormalization. 



and the absence of these counterterms uniquely fixes the coefficients in (|12.1|) . The explicit 
computation of the ZZ and WW wavefunction renormalizations demonstrates that the above 
operator has the correct relative factor between the two gauge bosons. This is perhaps non- 
trivial, since there is no additional degree of freedom to fix this relative factor. We also note 
that an identical operator to ( |12.2| ) is also required in technicolor theories in order to cancel 
the divergent contribution of pseudo-Goldstone bosons to the oblique parameters [^. 
The operator which provides the counterterms for the mass renormalization is 

Om = ^ (g'\D,H^H){H^D^H) + H {D H)) . (12.3) 

The first operator that appears here violates the custodial symmetry, and in particular 
contributes to the Z but not the W mass. After electroweak symmetry breaking they together 
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reduce to 




M^W;[W-^ + -^Z^Z^ . (12.4) 



Note that it is my which appears, so we exphcitly see that this operator contributes to 
the p parameter. In this case it is trivial to obtain the correct mass renormahzation from 
( |12.^ ), since here there are two coefficients to be determined from only two constraints. So 
in addition to the usual Standard Model renormalizations, these wavefunction and mass 
counterterms are also required to renormalize the model. 

The model with the radion represents an effective theory valid for E K. The dimension- 
6 operators discussed above are obtained by integrating some unknown degrees of freedom 
in the full theory, and in the effective theory they appear with some unknown coefficients 
Zj(A). These for example could include the effects of integrating out the heavy spin-2 KK 
modes. Since the divergences for which the two above operators act as counterterms arise 
at one loop order, they are proportional to 7^/(167r^). Thus it is reasonable to expect that 
the finite part of the operator is also of the same order, and in order to match the form of 
the explicitly calculated one-loop corrections we will write the (finite) coefficients as 

Zm — ^„ p p QM, Zx — p p Qx, (12.5) 

where we expect that the dimensionless parameters a^/ and ax are at most of order one. Note 
that since 7^ ~ A~^, in this parameterization the dimension-6 operators are still suppressed 
by A2. 

These coefficients parameterize the unknown physics integrated out at the scale ~ A. 
Since, however, the radion at one loop contributes to the anomalous dimension of these 
operators, when comparing to the experimental results evaluated at the Z mass large loga- 
rithms of 0{lnA/Mz) appear from this anomalous scaling and this effect should be included. 
Following Wilson, a one-loop Wilsonian renormahzation group equation is obtained for the 
operator coefficients. In the leading logarithm approximation the value of these coefficients 
at the weak scale is determined to be 

Z,(M^) = Z,(A) + -^ln^, (12.6) 

with the (3i determined from the coefficient of the In /i term (or equivalently, from the 1/e 
poles) in an explicit one-loop computation. To compute the oblique parameters, one then 
adds the contribution of these renormalized operators to the finite parts of the one-loop 
diagrams. In the leading logarithm approximation this amounts to simply replacing the 1/e 
poles in the gauge-boson self-energies with InA/M^. 

To compute the oblique parameters one uses the Feynman rules in the previous Section 
to compute the two Feynman diagrams that contribute to the vacuum polarizations. Both 
diagrams have one internal radion, and one uses the three point function and the other uses 
the four point function. Our convention for the sign of the vacuum polarizations Hvv is that 
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\ /wwvN^^^ywvwN Vv =in^V(P^) =^^''''nyy(p^)+2J9Vnyy(p2) (12.7) 



and only the first term is computed. The generic form of the radion contribution is 

Uvvip') = Ilvvip') + ^vvip')- (12.8) 

Here M' denotes the anomalous contribution due to the conformal coupling of the radion, 
and 'S" denotes the standard contribution which is also similar to the Higgs contribution 
(when ^ = 0). The anomalous couplings are discussed in the previous Section. By an 
appropriate rescaling of the coupling we can also use 11'^ for the Higgs. When ^ 7^ 0, the 
results given below can also be used to compute the oblique parameters after an appropriate 
redefinition of the couplings and masses. The modification to the expression for the oblique 
parameters is summarized in ( p.2.21| ). 

Inspecting the Feynman diagrams for the vacuum polarizations one finds that the quad- 
ratic divergences cancel between the two diagrams, leaving only a logarithmic divergence. 
Therefore this justifies the use of dimensional regularization. An explicit computation of the 
Feynman diagrams in unitary gauge and for vanishing curvature scalar parameter ^ gives 

s (6 5 ml m\\iim\ j 11^ - mlhimH 

iiyv-iuj - -T^— 1 7 + o - T^zri" + —2 — 1:2 ) ' l^^.yj 




mt — m 



V 



/ (ixln(x^my + (1 — x)m^)//i^ 

JO 



777^ 

.2.2 I ^2 .,^2 -l)ln-f + 7&-2)ln:^) , (12.10) 




IGvr^ ^2 3my 



where to avoid confusion with defining too many 7's, the renormalization scale ^ appearing 
includes the usual factors of 47r and Euler's constant 7^;. An analytic expression for the 
Feynman parameter integral may be obtained, but it is not very illuminating. A powerful 
check on these expressions is gauge invariance. We have explicitly checked that in the general 
gauge, the gauge parameter a cancels from the expression and reproduces the above 
results. We note that the divergences appearing here have the form as given by the operators 
in ( |12.1| ) and ( p.2.3| ). For H (0)^/1/ ~ my as required by ( |12.3| ). The difference between H(my) 
and H(0) gives the divergence proportional to p^, but in the above equations has already 
been set to my. With this in mind, by inspection the coefficient of is ~ my, as required 
by (|12.1|) . So the set of operators given by (|12.1| ) and (|12.3| ) provide an appropriate set 
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of counterterms. We renormalize using MS. Then using the Wilsonian approach outhned 
above, in the leading logarithm approximation we just replace 1/e with InA/M^, that is 



2 , . , 



^-7s + ln47r^ln^, (12.11) 



and set /x = Mz- 

The anomalous contribution is finite and is 



We note that the term does not contribute to the oblique parameters since it is p 
independent and also ~ m^. 



The PDG convention for the oblique parameters that we use here is 

_ 1 /iwo) n^z(o) \ 

^ 4sin^g^cos^g^ ( lizzjMl) n^z(0) \ . . 

a \ Ml Ml I ^ ■ > 

S + U = 4sin^^H/ ( UwwjM^) ^ww{0) \ .^2 15) 

a [ M^ M^ J ^ ■ ^ 

More generally one must also include the 11^^ and 11^^ self-energies. They have been dropped 
here since they do not receive contributions from either the Higgs or the radion. 

Using the above expressions one can evaluate the contribution of the radion to S and T 
in the limit of a large radion mass. One obtains for ^ = 0, 



TT V 12 mi 72 ' 



T = _i2;_(-^+,„4+l), (12.17) 

167rcos^6'vK \ 3 m; 6 y 

Inspecting the above expressions for the 11 's one finds that there is no divergent contribution 
to U, and this is consistent with the fact that Ox and Om provide no counterterms for U. 

It is interesting that the radion contribution to S is negative and to T is positive. This 
is easy to understand by comparing this result to the contribution of the Higgs in the SM. 
In fact, for 7 = 1, m,. = mh and Zi{A) = 0, the radion result is identical to the (logarithmic) 
contribution of the SM physical Higgs. But there the total correction to S and T is finite, so 
that for the Higgs the A dependence in the above expression is canceled by loops of H^'s and 
Z, leaving a In nih/Mz dependence. Then one obtains the usual positive (negative) correction 
to 5" (T). For the radion though the large logarithms are present, and since A > m^, the 
radion contribution to S is negative, and to T it is positive. Recalling that 7^ = v^/dA"^, the 
size of the above correction is only significant for small values of A. 
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We conclude with a comment on the decouphng behavior of the radion. Inspecting the 
above expressions we see that for large the radion contribution scales as 

1 Kr 

— In^. (12.18) 
h? rur 

For the purposes of this paragraph we distinguish the cutoff of the effective theory, Ap, 
from the mass scale A = Mpie~'^'^° appearing in the radion coupling. In our analysis we 
have approximated Ac ~ A. If the radion mass is much larger then A, then the cutoff of the 
effective theory is the much higher scale E ~ or larger. Then it is more natural to express 
all higher dimension operators as suppressed by m^^ or A^^. But then the coupling of the 
radion to the gauge bosons contains the very large coefficient m^/A. It is then inappropriate 
to use the one-loop approximation. So to remain within the validity of the approximations 
used here one must also increase A for large rrir- Then the radion contribution decouples. 

For large radion mass it is conceivable that the coupling of the radion to TrT decreases. 
This cannot be seen in our computations here because we have made the approximation 
of using only the zero mode wavefunction to determine the coupling, but have included 
the backreaction perturbatively to compute the radion mass. For very large radion mass 
these approximations are invalid, and then one must exactly solve the equations. It is then 
conceivable that for large radion mass the radion wavefunction on the TeV brane decreases 
in such a manner that the coupling to TrT remains natural, i.e. 0(m~^). 

Following standard practice we define a reference model in which one computes the 
oblique parameters within the Standard Model, which means for some specific value for 
the Higgs mass. Since the curvature scalar operator mixes the radion and Higgs, the two 
physical scalars are some mixture of the gauge Higgs and radion. This mixture is not a 
unitary rotation due to the kinetic mixing between the states. The couplings of the 'Higgs' 
in this case is somewhat different than in the Standard Model, and for the purposes of 
computing the oblique parameters it is easiest to think of this as a new model, rather than 
as a perturbation to the Standard Model. So in computing the oblique parameters the 
Standard Model Higgs contribution (for the reference Higgs mass) should be subtracted out, 
and the contribution of the physical states in this model added back in. That is, 

^ - ^SM = X^ewirrih: m^, A, ^) - X^^ {mu = rn^^). (12.19) 

As mentioned, X^^^rrih = rn^^) is the contribution from only the Higgs, with mass rn^^ and 
with Standard Model couplings, and it is independent of A, and the curvature-mixing 
parameter ^. The quantity X'^g^ is the full SM contribution, with the Higgs set at the 
reference mass rri^^ . The new physics contribution contains two pieces, 

^new = + , (12.20) 

which describe the contribution of the physical radion {Xf>), and the physical Higgs (Xh)- 
In the limit of no curvature-scalar mixing, this last contribution is just that of a Standard 
Model Higgs with mass ruh. For a general curvature scalar mixing one just needs to include 
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the effect of the mixing coefficients a, b, c and d given in the previous section. Then 
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Here X"^ is the anomalous contribution, obtained from the vacuum polarization (|12.12| ), but 
dropping the term since this does not contribute to any of the oblique parameters. The 
other X's are obtained from using the vacuum polarizations ( 12.9 ), ( 12.10 ), and inserting 
the physical mass of the appropriate state. 

We note that, for example, the full anomalous contribution to S from the above formula 

is 

' (6^-1) (12.22) 



A 

new 



QZ'^TT 

and is negligible for reasonable values of 7 and ^. The anomalous contribution to T is even 
smaller. 



12.1 Numerical Results 

The "new" contribution, X, is constrained to lie within the measured values (extracted 
assuming mf^ '""^ = 100 GeV) H 

^meas = -0.07 ± 0.11 (12.23) 

Trne.s = -0.10 ±0.14 (12.24) 
U^e^ = 0.11 ±0.15. (12.25) 

(The errors are for one sigma.) X can be easily calculated once rrih, rrir, A, and ^ are 
specified. Notice that since the current best-fit values for the electroweak parameters are 
nonzero, the new contributions can be more weakly or more strongly constrained depending 
on whether they add destructively or constructively with the Higgs, respectively. As a first 
example, we show the contribution to S and T in Fig. ^ as a function of ruh = rrir (the 
"gauge" masses), fixing A = 1 TeV. Each contour corresponds to a different value of ^, and 
the contours end when a physical mass exceeds the cutoff. The unshaded region corresponds 
to the la allowed region. Notice that T is a strong constraint on small (gauge) masses, while 
S* is a strong constraint for large masses. Also, the the case with ^ = is nearly identical to 
the ordinary SM Higgs contribution, since the radion contribution that can be separated out 
is strongly suppressed by the coupling 7^. In Fig. ^ we show the contribution to S and T as 
a function of rrir with rrih = 300 GeV. Notice that the contributions are nearly independent 
of rrir for small curvature scalar mixing. 

The above results illustrate a general trend that with small or absent curvature-scalar 
mixing, the bound on the Higgs mass is not significantly affected by the presence of the 
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Figure 3: The contributions to S*, T as a function of the "gauge" masses mh = rrij.. Each 
hne is a contour for a fixed curvature scalar mixing ^. The cutoff scale was chosen to be 1 
TeV (7 = 0.1). The shaded regions are excluded by the PDG measurements to one sigma. 

radion. This is not true, however, if we allow ^ to take larger values. It is easiest to first 
illustrate that radion physics with large curvature scalar mixing can significantly relax the 
bound on the Higgs mass, by scanning through the parameter space (choosing irih = rrir for 
simplicity) for values that satisfy one or two sigma limits on the electroweak parameters. We 
find sets of parameters that are not minor perturbations on the SM limit allow the physical 
masses of the Higgs and radion to be several hundred GeV, and perhaps even TeV scale. In 
Fig. 1^ we show the range of physical masses and the range of ^ as a function of the cutoff 
scale. In general there is not a unique mapping between the figures, however the "shark 
fin" structure for the one sigma region in Fig. ^(a) does correspond to the "inverse fin" in 
Fig. |^(b). Notice that at two sigma the physical Higgs mass can be much larger than the SM 
bound throughout the parameter space shown, and even at one sigma there exists a narrow 
range of large, negative curvature-scalar mixing where the physical Higgs mass could be of 
order a TeV. The latter result arises from a cancellation between the physical Higgs and 
radion contributions with the SM reference contribution. This can be seen in the limit of a 
large radion and Higgs mass. Since the dependence on the masses is only logarithmic, we 
can approximate the masses as being equal. Then, for example. 



The first contribution is just the usual correction from the Higgs. But the second correction 
can be potentially large and negative due to the dependence on the ^ parameter. It should 
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Figure 4: Same as Fig. ^ except that is fixed to 300 GeV. Notice that the contributions 
to S and T are nearly independent of the radion mass if the curvature mixing is small since 
the radion contribution is suppressed by ~ 7. 

be emphasized that the large correction is due to the (non-unitary) kinetic mixing between 
the radion and Higgs, or equivalently, due to the non-standard couplings of the radion and 
Higgs in the mass basis. Hence, while this region is provocative, it nonetheless requires 
fine-tuning. 

These results have assumed that the contribution from the nonrenormalizable countert- 
erms is small, meaning ax and are less than order 1. For larger coefficients the allowed 
regions of parameter space, albeit only at moderately large 7. In Fig. ^ we show the shift in 
the contours, for the two sigma region, resulting from taking ax = ±10. [a^ was also taken 
to be 10, but the effect on the contours was negligible.) 

13 Limits on Radion Mass 

As we found in Section H, the mass of the radion is expected to be significantly below the the 
cutoff scale, placing it in a region that can be directly probed by experiments. The previous 
section has shown that the radion couples much like a Higgs boson, and in the limit ^ 0, 
the tree-level couplings of the radion are simply scaled by 7. Let us first consider the bounds 
in this case. 

In the SM, the current bound on the Higgs mass comes primarily from the LEP processes 
e^e^ ^ Z* ^ Zh, with the value mf^ ^ 108 GeV |3^. For the radion, an exactly analogous 
production process occurs e^e^ ^ Z* ^ Zr, except that the ZZr coupling has a factor 
of 7. To a good approximation, we can therefore estimate the production cross section of 
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Figure 5: The allowed region of and as a function of the inverse of the cutoff scale 

7 = f /A by requiring S, T, U do not exceed the one sigma (dark region) or two sigma (light 
region) measurements from the PDG. The dashed lines correspond to the theoretical bound 
requiring the kinetic term is non-negative [see Eq. (9.14)]. The black sliver corresponds to 
the region where mf^^'^ ~ 300 GeV. 

radions at LEP by simply scaling the Higgs cross section by 7^. 

The decay of the radion is somewhat more complicated, however. As we discussed in 
Section [Tl|, the radion couples directly to gauge bosons through the conformal anomaly. 
Although this coupling is one-loop suppressed, it competes with Yukawa suppressed inter- 
actions and, for the case rgg, can be comparable or even dominate ||21|, In the radion 
mass range well below the tt threshold, the ratio of the two largest widths can be expressed 
as 



where = 1 — 4mf/m^ and C3 ~ 23/3 is roughly the one- loop QCD /3-function coefficient 
(approximately including the smaller contribution from the one-loop triangle diagram with 
top quarks.) Notice that the coupling 7^ cancels in this ratio. This can be written in the 
suggestive form 



Hence, r gg dominates for the region 12m;, ^ ^ 2Mw The search strategy for 
the radion is therefore significantly different from the Higgs in this mass window, namely 
searching for a pair of gluon jets instead of a pair of 6-jets. Similarly, the radion has a 
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Figure 6: The shift in the two sigma contours shown in Fig. § resulting from taking the 
coefficients of the nonrenormahzable operators to be ax = 10 (dark region) and ax = —10 
(hght region). 



different production cross section at hadron colhders via gluon fusion, proportional to the 
conformal anomaly-enhanced width into gluons but suppressed by the usual 7^ . 

Determining an accurate bound on the radion mass in the region that can be probed by 
LEP requires a detailed analysis of detecting a two gluon plus Z signal. We will not attempt 
this here. Instead, the expected bound on the radion mass can be roughly estimated as 
a function of the coupling if we assume that some number of production events at LEP 
could not have escaped detection (or be lumped into SM backgrounds) . Near the kinematical 
limit, the best bound will always come from the highest energy data. For lower mass radions, 
a lower center-of-mass energy results in a slightly higher cross section. [] Since the bound for 
a given radion mass is limited only by luminosity, we can combine the multitude of LEP 
runs at various energies by weighting by the integrated luminosity accumulated. The bound 
is then simply 

N 

(13.3) 



f < 



E crv^(e+e Zh; nih = rrir) x J C^, 

where the sum is over the various recent LEP runs with center-of-mass energy -^i. 
encodes all of the detailed analyses of backgrounds, signal efficiencies, etc., and is in general 
not independent of energy or radion mass. In Fig. |^ we simply show the bound obtained if 
= 20 or 100 (the integrated luminosity for each energy was also summed) corresponding 
to producing 20 or 100 events summed over all four LEP experiments. These numbers were 
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Figure 7: The bound on the radion mass as a function of 7 is shown, assuming the signal 
could be extracted from background once the radion production cross section times integrated 
luminosity exceeds = 20, 100 events at LEP (summed over the four experiments). 

chosen since searches for Higgs bosons typically need a few to tens of events (per detector) 
for a statistically significant signal-to-background ratio. Notice that no bound on the radion 
mass is expected from the recent LEP runs once 7 is less than about 0.1. 

One could also search for light radions, ^ 60 GeV, at LEP I via the decay Z —* ff + r, 
through the same coupling discussed above. However, the expected bound obtained by using 
this procedure is no better than that found above for larger than about 10 GeV. 

We have not attempted to estimate a bound on 7 for a radion mass less than about 10 
GeV. We really do not expect the radion to be several orders of magnitude below the cutoff 
scale, and so at the outset it seems this mass region is unnatural. But, the presence of several 
low energy production processes (and rare decays) could be important, so a considerably more 
careful analysis than what we have attempted here is needed. 

When curvature scalar mixing is included, the coupling ZZr is modified as shown in 
Eq. (|11.20|) . The above analysis can be translated into this more general case, but now the 



coupling is not simply 7 but a function of the curvature-scalar mixing as well. In addition, 
the SM Higgs couplings are also modified, and so its production and decay are also affected. 
In particular, the production cross section could be either enhanced or suppressed. (This is 
similar to what happens in two Higgs doublet models, such as the MSSM.) An interesting 
signal for RS with curvature-scalar mixing could be observing a nonstandard cross section 
or decay rate for a SM-like Higgs. 
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14 Conclusions 



In this paper we have analyzed the coupled radion-scalar system in detail, including the 
backreaction of the bulk stabilizing scalar on the metric. We derived the coupled differential 
equations governing the dynamics of the system, and found the mass eigenvalues for some 
limiting cases. We find that due to the coupling between the radion and the bulk scalar, 
there will be a single KK tower describing the system, with the metric perturbations non- 
vanishing for every KK mode. This implies that the Standard Model fields localized on the 
TeV brane will couple to every KK mode from the bulk scalar, and this could provide a 
means to directly probe the stabilizing physics. We also found that in an expanding universe 



the shift in the radion at late times completely agrees with the effective theory result of [|I4 

We also calculated the contributions of the radion to the oblique parameters using an 
effective theory approach. Since the radion is the only new state well below the TeV scale, we 
argued that a low-energy effective theory including only the radion and SM fields is sufficient, 
as long as appropriate nonrenormalizable counterterms at the cutoff scale are added. In the 
absence of a curvature-scalar Higgs mixing operator, the size of the contribution to the 
oblique parameters due to the radion is small. In the presence of such a mixing operator, 
the corrections can be much larger due to the modified radion and Higgs couplings. In 
particular, including only the mixed radion and Higgs fields as "new physics" , we calculated 
the range of curvature-scalar mixing for a given cutoff scale that allows the physical Higgs 
mass to be up to of order the cutoff scale, while S^ew and T^ew were within the experimental 
limits. However, the parameters must be increasingly fine-tuned to achieve a Higgs mass 
that exceeds a few hundred GeV. 
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